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Abstract 

In this paper, we investigate the dark matter halo density profile of M33. 
We find that the HI rotation curve of M33 is best described by a NFW dark 
matter halo density profile model, with a halo concentration of Cvir = 4.0 ±1.0 
and a virial mass of M^ir = (2.2 ± 0.1) x 10^^ M©. We go on to use the NFW 
concentration (cvir )of M33, along with the values derived for other galaxies 
(as found in the literature), to show that Cvir correlates with both spiral arm 
pitch angle and supermassive black hole mass. 
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1. Introduction 

The currently favored cosmological model, Lambda+ Cold Dark Matter 
(ACDM), is remarkably successful at reproducing the large-scale structure of 
the Universe (Blumenthal et al. 1984; Springel et al. 2005). However, small- 
scale observations have proven harder to explain. High-resolution N-body 
simulations of ACDM structure formation predict that the central density 
profiles of dark matter halos should rise steeply at small radii, p{r) cx r~'^, 
with 7 ^ 1 — 1.5 (Navarro, Frenk & White 1997, henceforth NFW; Navarro et 
al. 2004; Diemand et al. 2005). Observations of rotation curves of late- type 
disk galaxies and dwarf galaxies, on the other hand, have shown that quite 
often, mass distributions with lower than predicted densities or with constant 
density cores, where 7 ~ (i.e., a pseudo-isothermal profile), are preferred 
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(Swaters et al. 2003; Gentile et al. 2004, 2005; Simon et al. 2005; Kuzio de 
Naray et al. 2006, 2008; Shankar et al. 2006; Spano et al. 2008). This is 
known as the cusp/core problem. One possibility is that these observations 
are pointing to a real problem with ACDM cosmology, perhaps indicating 
that the dark matter is not cold, but rather warm (Zentner & Bullock 2002), 
in which case it is easier to produce constant density cores at the centers of 
dark matter halos. Another possibility is that these late-type galaxies have 
constant density cores because of their late formation (Wechsler et al. 2002) 
and that earlier-type bulge-dominated galaxies (which form at earlier times) 
will tend to conform to the standard expectations of the theory. This is 
because the central mass densities of galaxies tend to reflect the density of 
the Universe at their formation time (Wechsler et al. 2002). 

In this paper we have chosen to model the HI rotation curve of M33 
from Crobelli & Salucci (2000). Due to its proximity, M33 can be studied in 
exquisite detail, and it therefore provides a crucial testing ground of our ideas 
of galaxy formation. Its Hubble classification is SA(s)cd (de Vaucouleurs et 
al. 1991), meaning that is of particularly late-type, with little or no bulge. 
This is reflected in the central supermassive black hole mass of Mbh < 1500 
Mq (Gebhardt et al. 2001), and black hole masses tend to be related to the 
central bulge mass (Magorrian et al. 1998; Haring & Rix 2004). In this paper 
we model the rotation curve of M33 with both a pseudo-isothermal profile 
dark matter halo density model and an NFW dark matter halo density model. 
We then use parameters derived from these fits to look at relations between 
the dark matter halo and other galaxy properties, such as supermassive black 
hole mass and spiral arm pitch angle. 

This paper is organized as follows. Section 2 describes the observed data 
and data analysis. Section 3 describes how the rotation curve is modeled 
and how we derive the baryonic and dark matter halo contributions to the 
rotation curve. Section 4 discusses our results and Section 5 summarizes our 
findings. Throughout this paper, we assume a fiat ACDM cosmology with 
Qrn — 0.27 and a Hubble constant Hq — 75 km s~^ Mpc~^. 

2. Observations and Data Reduction 

We have made use of the Spitzer/IRAC 3.6-//m image of M33. The IRAC 

observations were taken as part of the Gehrz Guaranteed Time Observer 
Program ID 5. The mapping sequence for each epoch consisted of ^ 148 
positions per channel. Each position was observed with three 12 s frames 
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dithered with the standard, small, cycling pattern. The FWHM of the point- 
spread function (PSF) at 3.6-/im is 1.7" or 6.9 pc at the distance of M33. 
The final mosaic spans an area of ~ l.'O x 1.'2. We adopt a distance to M33 
of (i = 840 kpc (e.g., Magrini, Corbelli & Galh 2007), and it has a redshift 
of z = -0.000597 (de Vaucouleurs et al. 1991). 

For dynamical measurements, we make use of the HI rotation curve of 
Corbelli & Salucci (2000). We also make use of the inclination corrected HI 
linewidth from HyperLeda0 of 100.4 ± 3.0 km s~^ (e.g., Paturel et al. 2003). 

For the determination of the spiral arm morphology we have made use of 
an R band image from the Digital Sky Survey (DSS). 

2.1. Measurement of spiral arm pitch angle 

Spiral arm pitch angles are measured using a two-dimensional fast Fourier 
decomposition technique, which employs a program described in Schroder et 
al. (1994). Logarithmic spirals are assumed in the decomposition. 

The amplitude of each Fourier component is given by 

ELi E/=i ^ii (In r,e)exp- [2 (m^p In r)] 

^{^,P) - Tn Q\ 

Ei=i Ej=i^ij(lnr,^) 

where r and 6 are polar coordinates, J(lnr, 6') is the intensity at position 
(Inr, ^), m represents the number of arms or modes, and p is the variable 
associated with the pitch angle P, defined by P = —{m/pmax)- Throughout 
this work we measure the pitch angle P of the m = 2 component. 

Pitch angles are determined from peaks in the Fourier spectra, as this is 
the most powerful method to find periodicity in a distribution (Considere & 
Athanassoula 1988; Garcia-Gomez & Athanassoula 1993). 

The image was first projected to face-on. Mean uncertainties of position 
angle and inclination as a function of inclination were discussed by Considere 
& Athanassoula (1988). For a galaxy with low inclination, there are clearly 
greater uncertainties in assigning both a position angle and an accurate in- 
clination. These uncertainties are discussed by Block et al. (1999) and Seigar 
et al. (2005, 2006), who took a galaxy with low inclination (< 30°) and one 
with high inclination (> 60°) and varied the inclination angle used in the 
correction to face-on. They found that for the galaxy with low inclination, 
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the measured pitch angle remained the same. M33 has a relatively low in- 
clination of ~ 30°, and so the uncertainty in the inclination angle in this 
case, does not result in a large error in the pitch angle we measure for M33. 
Our deprojection method assumes that spiral galaxy disks are intrinsically 
circular in nature. 

3. Mass modeling 

3.1. The baryonic contribution 

Our goal is to determine a mass model for M33 from direct fitting of mass 
models to its rotation curve. We perform a bulge-disk decomposition in order 
to estimate the baryonic contribution. We then determine several different 
models and try to recreate the nuclear spiral by minimizing reduced-x^- 

We first extract the surface brightness of M33 using the Spitzer 3.6-fim 
image and the IRAF ELLIPSE routine, which fits ellipses to an image using 
an iterative method described by Jedrzejewski (1987). In order to mask 
out foreground stars, SEXTRACTOR (Bertin & Arnouts 1996) was used. An 
inclination correction was then applied to the surface brightness profile (de 
Jong 1996; Seigar & James 1998) as follows 



where /x, is the surface brightness when viewed at some inclination z, /i is the 
corrected surface brightness, a is the major axis, b is the minor axis and C is 
a factor dependent on whether the galaxy is optically thick or thin; if C = 1 
then the galaxy is optically thin; if C = then the galaxy is optically thick 
(e.g., Seigar & James 1998; de Jong 1996). Graham (2001a) showed that 
C = 0.91 is a good value to use for the near-infrared Kg band. Adopting 
a simple reddening law, where extinction falls as the square of wavelength, 
it can be shown that a value of C = 0.97 is appropriate at 3.6-/im (Seigar, 
Barth & Bullock 2008a) and we adopt this value here. 

The resulting surface brighntess profile Fig. [T] reaches a surface bright- 
ness of /i3.6 ~ 20.7 mag arcsec"^ at a radius of ~13.2 kpc (equivalent to 54.0 
arcmin). From this surface brightness profile, we perform a one-dimensional 
bulge-disk decomposition, which employs the Sersic model for the bulge com- 
ponent and an exponential law for the disk component (e.g., Andredakis, 
Peletier Sz Balcells 1995; Seigar & James 1998; Khosroshahi, Wadadekar & 
Kembhavi 2000; D'Onofrio 2001; Graham 2001b; Mdllenhoff & Heidt 2001; 
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, fj,^ = 19.57 mag arcsec"' 

E = 0.39 kpc 

1' - n = 1.0 
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Figure 1: The Spitzer 3.6-/im surface brightness profile with decomposition into bulge and 
disk components. The bulge has been fitted with a Sersic model (short-dashed line) and 
the disk has been fitted with an exponential model (long-dashed line). 



Table 1: M33 Observational data. The Hubble type is from de Vaucouleurs et al. (1991). 
The distance in kpc is taken from Magrini et al. (2007). 



Parameter 



Measurement 



Hubble Type 
Distance (kpc) 

Position angle of major axis (°) 
Bulge effective radius, Re (arcmin) 
Bulge effective radius, Re (kpc) 

Bulge surface brightness at the effective radius, fXe (3.6 /xm-mag arcsec"^) 
Bulge Sersic index, n 

Disk central surface brightness, iJ.o (3.6 /zm-mag arcsec"^) 
Disk scalelength, h (arcmin) 
Disk scalelength, h (kpc) 
Disk luminosity, Ldisk (L©) 
Bulge-to-disk ratio, B/D 



SA(s)cd 

840 

23 

1.60±0.11 
0.39±0.03 
19.57±0.98 
1.0 

18.08±1.02 

6.95±0.49 

1.70±0.12 

(3.16 ± 0.30) X 10^ 

0.03 
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see Graham & Driver 2005 for a review). The Sersic (1963, 1968) i?^/'^ model 
is most commonly expressed as a surface brightness profile, such that 



Ijl{R) = He exp I -bn 



Re) 



(3) 



where fi^ is the surface brightness at the effective radius that encloses half 
of the total light from the model (Ciotti 1991; Caon, Capaccioli & D'Onofrio 
1993). The constant 6„ is defined in terms of the parameter n, which de- 
scribed the overall shape of the light profile. When n — A, the Sersic model 
is equivalent to a de Vaucouleurs (1948, 1959) i?^/^ model and when n = 1 it 
is equivalent to an exponential model. The parameter bn has been approxi- 
mated by bn = 1.9992n- 0.3271, for 0.5 < n < 10 (Capaccioli 1989; Prugniel 
& Simien 1997). The exponential model for the disk surface brightness profile 
can be written as follows 



A((i?) = //oexp(-i?//i) (4) 

where /iq is the disk central surface brightness and h is the disk exponential 
scalclcngth. The results of our surface brightness fitting are summarized in 
Tabic 1. 

We now assign masses to the disk and bulge of M33. The stellar mass- 
to-light ratio in the Kg band is a well-calibrated quantity (Bell et al. 2003) 
which depends on B — R color. Seigar et al. (2008a) extended this to a 3.6- 
/xm image of M31 using the population synthesis codes of Bruzual & Chariot 
(2003) and Maraston (2005). Using their results, we find a central mass-to- 
light ratio of M/Lg.g ~ 1.25 ± 0.10 with a gradient of -0.014 kpc"^ This 
results in a disk mass of Mdisk = (3.81 ± 0.47) x 10^ and a bulge mass of 
Mbuige = (1.14 ± 0.14) X 10^ M© for M33. 

A concern in using the 3.6-/xm S'/Jzteer waveband to determine the under- 
lying stellar mass, is the effect of emission from hot dust in this waveband, 
although this is probably only important in or near Hll regions. In order to 
place some constraint on this, we have chosen to explore the emission from 
dust in the near-infrared K band at 2.2 //m. Using near-infrared spectroscopy 
at 2.2 //m, it has been shown that hot dust can account for up to 30 per cent 
of the continuum light observed at this wavelength in areas of active star 
formation, i.e., spiral arms (James & Seigar 1999). When averaged over the 
entire disk of a galaxy, this reduces to a 2 percent effect, if one assumes that 
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Figure 2: The HI rotation curve from Corbelli & Salucci (2000) modeled using a pseudo- 
isothermal model (core model; blue solid line) and a NFW model (red dotted line). The 
squares represent the total rotation velocities, and the circles represent the contribution of 
the dark matter to the rotation velocities (after subtraction of the stellar and gas mass). 



spiral arms can be up to 12° in width. At 3.6 /xm, this would therefore result 
in 3 percent of emitted light from dust. 

Another concern for the 3.6-/xm waveband would be the contribution from 
the polycyclic aromatic hydrocarbon (PAH) emission feature at 3.3 fim. How- 
ever, an Infrared Space Observatory (ISO) spectroscopic survey of actively 
star-forming galaxies by Helou et al. (2000) found that the 3.3-/im feature 
was very weak when they analysed the average 2.5-1 1.6-/im spectrum of 45 
galaxies. The contribution of the PAH feature to the 3.6-/im Spitzer wave- 
band, is therefore not a major concern. 

One other important contribution to the baryonic mass of M33 is the 
gas mass. Corbelli & Salucci (2000) have shown that beyond a radius of 10 
kpc, the gas contributes about the same to the rotation curve as the stars. 
Since, the best current estimate of the gas distribution comes from Corbelli 
& Salucci (2000), we have chosen to adopt their model for the distribution 
of gas mass in M33. 

3.2. The dark halo contribution 

A range of allowed dark matter halo masses and density profiles is now 
explored, using two models for dark matter halo density profiles, the pseudo- 
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isothermal model (e.g., Simon et al. 2005; see Equation and the Navarro, 
Frenk & White (1997; hereafter NFW) profile. A pseudo-isothermal density 
profile is given by 

which in terms of rotational velocity becomes 

K^(i?) = K^(oo)(l-§tan-i^), (6) 

where Rc is the core radius, and po = V^{oo) /Ai^GFi^. The NFW profile is 
given by 

P^^^ = {r/r,){iIr/r,Y 

where Rs is a characteristic 'inner' radius, and p° is the present critical density 
and 5c a characteristic overdensity. This overdensity is defined as 

lOOc^- 

5c = (8) 
where Cvh=Rvir/Rs is the concentration parameter and 

In (1 + Cvir) - Cvir/ (1 + C^ir) 

The circular velocity associated with this density is given by Battaglia et al. 
(2005) and is 



ta2 _ ^vhOjCy 
' c 



In (1 + CvirS) 



1 + Cv 



(10) 



where Kir is the circular velocity at the virial radius i?vir and s = R/R^ij-. 
This NFW profile is a two parameter function and completely specified by 
choosing two independent parameters, e.g., the virial mass M^ir (or virial 
radius -Rvir) and concentration Cyir = R^ir/Rs (see Bullock et al. 2001a for a 
discussion). Similarly, given a virial mass Mvir and the dark matter circular 
velocity at any radius, the halo concentration Cvir is completely determined. 

We now proceed by finding the best-fitting NFW and pseudo-isothermal 
(or constant density core) dark matter halo density profiles that describe the 
complete HI rotation curve of M33 as observed by Corbelli & Salucci (2000). 
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Table 2: M33 rotation curve modeling results, showing the best-fitting NFW and pseudo- 
isothermal models. 



Parameter NFW model 

'M~ (2.2 ±0.1) X 10^1 M, 

Cvir 4.0 ±1.0 

xVf^ 1-18 

Parameter Core model 

'R'c 1.39 ± 0.04 kpc 

1/(00) 105.4 ± 6.1 km s-^ 

XV/^ 3.19 



The result of this is shown in Figure 2. The pseudo-isothermal fit is shown 
as the solid blue line, with best-fitting parameters of ^(oo) = 105.4 ±6.1 
km and Rc — 1.39 ± 0.04 kpc, and a reduced-x^ value of x^/a* = 3.19, 
where ^ is the degrees of freedom. The NFW fit is shown as a dotted red 
line, with best-fitting parameters Cvir = 4.0±1.0 and Mvir = (2.2±0.1) x 10^^ 
M0, with a reduced-x^ value of x^/t^ = 1.18. As can be seen from Figure 
2, the pseudo-isothermal model (or core model in the figure) underestimates 
the rotation velocities beyond ~7 kpc. However, the NFW fit more closely 
recreates the observed data. This is also clear from the values of reduced- 
X^. We therefore conclude that the NFW model best represents these data, 
and this is consistent with the results of Corbelli & Salucci (2000). This is 
somewhat surprising for a late-type, bulgclcss galaxy like M33, since these 
late- type galaxies are often shown to have constant density cores (e.g., Kuzio 
de Naray, 2006, 2008). 

Table 2 lists the best-fit parameters of the best-fit NSF and pseudo- 
isothermal models based upon direct fitting to the HI rotation curve data. 

It is probably worthwhile noting that our best-fitting NFW model yields 
a concentration parameter, Cyir = 4.0 ± 1.0. This is somewhat lower than the 
concentration parameter of Cyir = 5.6 reported by Corbelli & Salucci (2000). 
Furthermore, we derive a virial mass of M^[y. = (2.2 ±0.1) x 10^^ M©, which 
is significantly lower than the virial mass of Mvu- = 7.4 x 10^-*^ Mq found by 
Corbelli & Salucci (2000). Here we discuss some reasons that could account 
for these apparent differences. Since we use the same gas distribution as 
Corbelli & Salucci (2000), the only difference can come from the stellar mass 
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component. The main difference l^etween our stellar mass component, and 
that of Corbelh & Salucci (2000), is that ours is determined from a Spitzer 
3.6-/xm observed in 2007, and that of Corbelh & Salucci (2000) is determined 
from a K band image reported by (Regan & Vogel 1994). The K band image 
from 1994 was taken when near-infrared arrays were really in their infancy, 
and so it is probably more important to rely on the more modern datasets 
when possible. Furthermore, Corbelh & Salucci (2000) assume a distance to 
M33 of 0.7 Mpc, whereas we use the more accurate measurement of 0.84 Mpc 
from Magrini et al. (2007). As a result of this underestimate in the distance 
to M33, Corbelli & Salucci (2000) have underestimated the size of the visible 
galaxy by a factor of ~17 percent, and this in turn has probably affected 
the total mass of M33 that they derive. Taking into account the different 
distances to M33, the disk scalelength of = 1.2 kpc used by Corbelli & 
Salucci (2000) would become h = 1.4 kpc if they had used the more accurate 
distance of 0.84 Mpc. This is still lower than the scalelngth oi h — 1.7 kpc 
that we report here. In converting this light distribution into stellar mass, we 
have then used a combination of the stellar mass-to-light ratios from Bell et 
al. (2003) and the population synthesis codes from Maraston (2005). These 
papers provide the best estimates currently available for determining the 
stellar mass-to-light ratios, and they were not available to Corbelli & Salucci 
when they performed their analysis. One final difference between our results, 
and those of Corbelli & Salucci (2000), is that we include the bulge mass, 
although considering the bulge-to-disk ratio of B/D = 0.03 this is unlikely 
to have a significant effect on the mass models. As a result, we conclude that 
the differences between our results and those of Corbelli & Salucci (2000), 
are caused by the different treatment of the disk starhght, updated stellar 
mass-to-light ratios, and more recent data. 

Finally, it should be noted that Corbelli & Walterbos (2007) revealed 
that M33 has a weak central bar. This could potentially have the affect 
of inducing non-circular motions in the central regions, i.e., within 1 kpc. 
However, Kuzio de Naray & Kaufmann (2010) have shown that, even in the 
case of barred galaxies, it is difficult to confuse an NFW dark matter halo 
profile with that of a pseudo-isothermal profile. In other words, our result 
that M33 is best described by an NFW profile, still holds, and given that the 
potential of the stellar bar is weak, the concentration is unhkely to change 
significantly. 

In the following discussion, we use the NFW concentration parameter to 
reveal some interesting relationships. 
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Figure 3; Spiral arm pitch angle versus rotation curve shear, showing a strong correlation. 
The solid squares represent galaxies with data measured by Block et al. (1999), the red 
squares are galaxies from Seigar et al. (2005), the blue squares are galaxies from Seigar 
et al. (2006), the cyan square is for Malin 1 (Seigar 2008), the green square is for M31 
(Seigar et al. 2008a) and the magenta square represents the data for M33 (this paper). 



4. Discussion 

Seigar et al. (2004, 2005, 2006) have demonstrated that a relationship 
exists between spiral arm pitch angle and rotation curve shear. Rotation 
curve shear is defined as: 

CO 2\ VdRj ' ^ ^ 

where A is the first Oort constant, u is the angular velocity, and V is the 
velocity measured at radius R. Using this equation it is possible to determine 
the shear from a rotation curve. We have performed such an analysis on the 
HI rotation curve for M33 and found a value for its shear of S = 0.46 ±0.01. 
We have also measured the spiral arm pitch angle for M33, which turns out 
to be P = 42.°2 ± 0.°3 (Seigar et al. 2008b). This pitch angle is in good 
agreement with previous measurements (Sandage & Humphreys 1980; Block 
et al. 2004). Figure 3 shows the relationship between spiral arm pitch angle 
and rotation curve shear. One can easily see that the pitch angle and shear 
values for M33 are consistent with the overall relationship. 
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Table 3: Spiral arm pitch angles, NFW concentration parameters, and central supermas- 
sive black hole mass for 5 galaxies. For the Malin 1 the spiral arm pitch angle; is taken 
from Moore & Parker (2006). For M31 the pitch angle is the average of values taken from 
Arp (1964) and Braun (1991). The NFW concentration value is taken from (1) Seigar 
(2008), (2) Klypin et al. (2002), (3) Seigar ct al. (2008a), (4) Seigar et al. (2006). The 
black hole mass estimates are taken from (5) Ghez et al. (2005), (6) Bender et al. (2005), 
(7) Gebhardt et al. (2001). 



Galaxy name 


Spiral arm pitch angle 
(degrees) 


Cvir 


Source 


(Mo) 


Source 


Mahn 1 


25.0 ± 1.0 


8.0 ± 1.0 


(1) 






Milky Way 




12.0 


(2) 


(3.7 ±0.2) X 10^ 


(5) 


M31 


7.1 ±0.4 


20.0 ±1.1 


(3) 


(1.7 ±0.6) X 10^ 


(6) 


M33 


42.2 ± 3.0 


4.0 




< 1500 


(7) 


IC2522 


38.8 ± 1.6 


8.0 ± 1.0 


(4) 






ES0582G12 


22.6 ±0.6 


22.0 ±5.0 


(4) 







Given the spiral arm pitch angles of a number of other galaxies, we can 
also now compare this quantity with the NFW concentration parameters for 
the galaxies hsted in Table 3. Figure 4 shows a plot of NFW concentration 
as a function of spiral arm pitch angle in degrees. This plot may only be 
for 5 galaxies, hut a relatively strong correlation appears to exist between 
these two quantities. Indeed Pearson's linear correlation coefficient is 0.95 
for this plot, although the significance at which the null hypothesis of zero 
correlation is disproved in onlt 54 percent, probably due to low number statis- 
tics. Nevertheless, an interesting correlation seems to exist between spiral 
arm morphology and dark matter concentration, and this could be further 
studied by targeting more galaxies in an observational campaign. Indeed, 
these data seem consistent with the suggestion that pitch angle and mass 
concentration are related (Seigar et al. 2005, 2006). 

Finally Figure 5 shows a plot of supermassive black hole mass as a func- 
tion of NFW concentration parameter. Unfortunately, here we only have 
data for three galaxies. Nevertheless, a hint of a correlation is starting to 
show, and seeing that such a correlation has been suggested by Seigar et al. 
(2008b), as well as Satyapal et al. (2008) and Booth & Schaye (2010), this 
plot is somewhat intriguing. This hint of a correlation should, of course, be 
expanded on by studying more galaxies along the Hubble sequence from type 
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Figure 4: NFW concentration parameter versus spiral arm pitch angle, showing a correla- 
tion. The green point represents data for M31, the cyan point for Malin 1, the red point 
for IC2522, the blue point for ES0582G12 and the magenta point shows the data for M33. 

Sa to Sd. 

5. Summary 

We have shown that the HI rotation curve of M33 can be best modeled 
with a dark matter halo that follows a NFW profile, with low NFW con- 
centration of Cvir = 4.0. Using the NFW concentration parameter from this 
fit, we find that interesting correlations between (1) spiral arm pitch angle 
and NFW concentration and (2) central supermassive black hole mass and 
NFW concentration, start to appear. Although the second correlation is only 
for three galaxies, on the surface it appears to be in disagreement with the 
argument made by Kormendy & Bender (2011) that the dark matter halos 
of galaxies have no affect on the masses of supermassive black holes found in 
their centers. These correlations are very intriguing and our results warrant 
further investigation, as we have been limited to data that was available for 
just a few galaxies. 
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